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Supercritical  water  oxidation  (SCWO)  in  hydrothermal  flame  regime  has  advantages  over  the  oxida¬ 
tion  in  flameless  regime.  The  main  advantage  is  that  the  feed  can  be  injected  into  the  reactor  at  low 
temperatures,  avoiding  plugging  and  corrosion  problems  in  a  preheating  system.  However,  there  is 
a  lack  of  kinetic  data  capable  of  properly  describing  the  flame  regime  oxidation.  In  this  study,  new 
global  reaction  rate  parameters  for  the  oxidation  of  isopropyl  alcohol  in  hydrothermal  regime  were 
adjusted  from  temperature  profiles  of  our  group’s  previous  experimental  data.  The  kinetics  obey  first 
order  rate  with  regard  to  the  fuel  and  the  oxidant,  and  follows  the  Arrhenius  law.  The  parameters  are 
ko  =  (9.308  ±  3.989)  x  1 07  (m3  s_1  kmol-1 )  and  Ea  =  89.441  ±  2.457  (kj  mol-1 ),  and  the  least  square  error  of 
the  fitting  was  1 0.8%.  This  kinetic  model  was  applied  in  a  parametric  analysis  of  flame  formation,  and  it  was 
used  to  analyze  the  behavior  of  a  supercritical  water  oxidation  vessel  reactors.  The  kinetic  model  is  able 
to  describe  the  behavior  of  the  vessel  reactor  when  working  in  steady  state  hydrothermal  flame  regime  at 
subcritical  injection  temperatures  (246  °C).  The  model  predicts  both  flameless  and  hydrothermal  flame 
regimes. 

©  2013  Elsevier  B.V.  All  rights  reserved. 


1.  Introduction 

Supercritical  water  oxidation  (SCWO)  in  hydrothermal  flame 
regime  is  a  promising  technology  for  the  total  destruction  of  waste 
because  overcomes  the  challenges  that  have  delayed  the  successful 
and  profitable  commercialization  of  this  technology:  corrosion,  salt 
deposition  and  high  energy  demand  [1,2].  Working  with  residence 
times  below  1  s,  the  process  can  be  performed  in  microreactors. 
When  the  temperature  of  the  mixture  becomes  greater  than  the 
ignition  temperature  a  flame  is  produced,  and  in  this  kind  of 
dense  aqueous  environment  it  is  known  as  hydrothermal  flame 
[3].  Hydrothermal  flames  were  referred  to  for  first  time  by  Franck 
and  coworkers  [4].  Several  research  groups  have  developed  reac¬ 
tors  working  with  a  hydrothermal  flame  as  heat  source  [5,6],  since 
SCWO  reactors  using  a  hydrothermal  flame  have  a  number  of 
advantages  [3]: 

•  Pollutants  can  be  destroyed  in  residence  times  of  a  few  milli¬ 
seconds. 

•  The  higher  operation  temperatures  improve  energy  recovery. 
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•  It  is  possible  to  initiate  the  reaction  with  feed  injection  tempera¬ 
tures  near  to  room  temperature  [7]. 

This  last  point  presumes  an  advantage  from  the  operational  and 
energetic  integration  points  of  view,  as  it  avoids  problems  of  plug¬ 
ging  and  corrosion  in  a  preheating  system. 

The  High  Pressure  Process  Group  (HPPG)  of  the  University  of 
Valladolid  (UVa)  (Spain)  recently  demonstrated  the  formation  of 
hydrothermal  flames  in  tubular  reactors  [8].  Using  this  simple 
device,  the  extinction  temperatures  of  the  flame  could  not  be 
reduced  to  less  than  370  °C.  Vessel  reactors  have  shown  to  be  more 
successful  in  maintaining  steady  stable  hydrothermal  flames  with 
injection  temperatures  near  to  room  temperature  [5,6,9].  The  flame 
stability  is  related  to  the  injection  temperature  and  flow  velocity, 
and  vessel  reactors  provide  flow  velocities  that  are  compatible  with 
hydrothermal  flame  front  velocities  [10].  It  is  thought  that  vessel 
reactors  provide  a  space  for  recirculation  in  which  the  cold  reagents 
are  preheated  to  the  ignition  temperature  and  are  brought  into  con¬ 
tact  with  the  radicals  already  formed,  making  it  possible  for  the 
flame  to  be  formed.  CFD  simulations  of  these  reactors  show  the 
existence  of  a  flame  front  where  most  of  the  reaction  takes  place 
[11]. 

This  development  on  the  SCWO  process  has  been  possible,  in 
part,  as  a  result  of  the  increase  in  the  quantitative  understanding 
of  the  SCWO  reaction  initiation  rates  of  different  types  of  organic 
species  at  the  lower  temperatures  of  450-500  °C  [12].  The  focus  of 
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Nomenclature 

C  molar  concentration  (kmol  m-3 ) 

Ea  activation  energy  (kj  mol-1 ) 

A Hr  heat  of  reaction  (J  kmol-1 ) 

M  molecular  weight  (kg  kmol-1) 

R  universal  gas  constant  (J  mol-1  K-1 ) 

Rk  species  source  (kmol  m-3  s-1 ) 

S  cross  section  (m2) 

Sh  energy  source  (J  m-3  s-1 ) 

T  temperature  (K) 

V  molar  volume  (m3  mol-1 ) 

X  auxiliary  parameter 

cp  constant  pressure  heat  capacity  (J  kg-1 1<-1 ) 

f0bj  objective  function 

h  specific  enthalpy  (J  kg-1 ) 

/<o  Arrhenius  pre-exponential  factor  (s  m3  kmol-1 ) 

fn  mass  flow  (kgs-1) 

r  reaction  rate  (kmol  m-3  s-1 ) 

t  time  (s) 

u  velocity  (ms-1) 

vt  translated  volume  (cm3  mol-1 ) 

w  mass  fraction 

y+  dimensionless  wall  distance 

z  longitudinal  coordinate  (m) 

Greek  symbols 

X  thermal  conductivity  (W  m-1 1<-1 ) 

p  density  (kg  m-3) 

r  stress  tensor  (Nm-2) 

Subscripts 
itjtl  directions 

k  species 

t  turbulent 

Abbreviations 

CFD  computational  fluid  dynamics 

Cl  confidence  interval 

PR  Peng-Robinson  equation  of  state 

IPA  isopropyl-alcohol 

SCW  supercritical  water 

SCWO  supercritical  water  oxidation 

TOC  total  organic  carbon 

TWR  transpiring  wall  reactor 

VTPR  volume  translated  Peng-Robinson  equation  of  state 


these  kinetic  studies  has  typically  been  on  model  compounds  rather 
than  on  actual  waste.  These  compounds  have  been  chosen  because 
they  contain  key  functional  groups  of  importance  and/or  repre¬ 
sent  the  rate  limiting  step  in  the  breakdown  of  a  range  of  complex 
waste  species.  Considering  only  homogeneous  reactions  without 
catalysts,  extensive  data  on  SCWO  kinetics  have  been  reported  for 
a  number  of  model  compounds  [12-16].  However,  these  kinetic 
data  from  literature  are  not  able  to  describe  flame  regime  oxidation 
[8].  Concentration  profiles  in  the  SCWO  reactor  are  very  difficult  to 
obtain  due  to  fast  reaction  rates  and  usually  only  global  conversions 
are  reported,  except  when  optical  devices  are  used  [17-19].  The 
lack  of  kinetic  models  is  one  of  the  reasons  why  most  of  the  model¬ 
ing  of  hydrothermal  flames  is  based  on  mixing  controlled  reactions 
[11,20].  This  approach  is  valid  for  non-premixed  flames  at  high  tem¬ 
peratures  since  it  assumes  that  chemical  reaction  is  fast  compared 
to  the  transport  processes  in  the  flow.  The  characteristics  of  our 


system  are  similar  to  premixed  flames,  where  fuel  and  oxidant  are 
totally  mixed  and  chemical  reaction  is  important. 

In  this  work,  a  new  kinetic  model  for  the  supercritical  water  oxi¬ 
dation  of  isopropyl  alcohol  in  flame  regime  has  been  adjusted  using 
our  group’s  previous  experimental  data  [8].  Temperature  profiles 
were  used  to  fit  kinetic  parameters,  by  developing  a  mathemati¬ 
cal  model.  This  model  is  able  to  reproduce  experimental  data  using 
either  air  or  oxygen  as  the  oxidant.  A  parametric  analysis  showed 
that  the  model  is  in  accordance  with  the  work  of  Serikawa  [21]. 
The  adjusted  kinetic  model  has  been  used  to  make  an  analysis  of 
the  behavior  of  a  vessel  reactor,  in  particular  of  the  transpiring 
wall  reactor  with  an  internal  hydrothermal  flame  studied  at  the 
University  of  Valladolid. 


2.  Experiments 

Experimental  data  used  for  adjusting  the  kinetic  parameters  of 
the  oxidation  of  isopropyl-alcohol  at  hydrothermal  flame  regime 
were  performed  at  the  University  of  Valladolid’s  SCWO  pilot  plant, 
which  had  been  adapted  to  work  with  tubular  mixers  as  described 
by  Bermejo  et  al.  [8].  The  pilot  plant  has  a  maximum  treatment 
capacity  of  20  kg/h  feed  and  uses  air  as  an  oxidant  compressed  by 
a  four  staged  reciprocating  compressor.  Both,  air  and  feed  are  elec¬ 
trically  preheated  before  being  introduced  into  the  tubular  reactor, 
which  consists  of  a  commercial  Ni-alloy  tube  of  1/4"  diameter  and 
1.5  m  length.  The  reaction  was  monitored  by  measuring  tempera¬ 
tures  at  the  mixing  point  and  at  several  points  inside  the  reactor, 
and  also  by  determining  the  total  organic  carbon  (TOC)  of  the  efflu¬ 
ent.  After  leaving  the  reactor  the  reaction  mixture  is  quenched.  The 
products  are  cooled  down  to  room  temperature  in  the  intercoolers, 
and  after  depressurization,  samples  of  the  liquid  and  gas  effluents 
can  be  taken. 

Additional  experiments  using  pure  02  as  the  oxidant  were  per¬ 
formed  in  a  similar  tubular  reactor,  installed  in  a  demonstration 
plant  of  the  University  of  Valladolid  situated  in  the  premises  of 
the  CETRANSA  company  in  Santovenia  de  Pisuerga  (Valladolid). 
In  this  case,  cryogenic  02  was  pumped  until  reaching  operational 
pressure  and  brought  to  room  temperature  by  an  evaporator.  The 
reactor  was  a  Ni-alloy  tube  of  1/4"  diameter  and  2  m  length.  The 
demonstration  plant  is  already  described  in  previous  work  [22]. 

To  determine  the  kinetic  parameters,  data  from  24  experiments 
using  air  and  8  using  02,  under  different  inlet  conditions  were  used. 

Once  the  kinetic  model  was  obtained,  it  was  introduced  in  a 
previous  model  of  the  transpiring  wall  reactor  working  at  the  same 
University  of  Valladolid  pilot  plant  as  described  by  Bermejo  et  al. 
[23],  and  the  results  of  simulations  were  compared  to  experimen¬ 
tal  data.  The  transpiring  wall  reactor  consists  of  a  stainless  steel 
pressure  shell  with  a  volume  of  10  L.  It  contains  a  reaction  chamber 
surrounded  by  a  porous  wall  through  which  clean  water  circulates. 
The  feed  and  the  air  are  introduced  into  the  reactor  through  its 
lower  part,  and  they  are  fed  through  the  injector  into  the  upper 
part  of  the  reaction  chamber;  the  reagents  flow  downwards  mixing 
with  the  clean  water  that  enters  the  reactor  through  the  transpir¬ 
ing  wall,  and  decontaminated  water  leaves  the  reactor  through  its 
lower  part.  More  constructive  details  of  the  reactor  can  be  found 
elsewhere  [10,23].  In  order  to  follow  the  reaction,  the  temperature 
is  measured  at  several  points  in  the  reaction  chamber. 


3.  Modeling 

3.2.  Tubular  reactor 

A  simple  model  of  the  tubular  reactor  described  in  Section  2  was 
programmed  in  order  to  fit  the  kinetic  parameters  of  IPA  oxidation, 
by  minimizing  the  difference  between  experimental  and  calculated 
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Table  1 

Volume  translation  at  23  MPa  for  application  in  VTPR-EoS. 
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vt  (cm3  mol-1 ) 

h2o 

4.291073 

o2 

-2.545861 

n2 

-3.743775 

co2 

1.284445 

IPA 

3.721115 

Fig.  1.  3D  tee-junction  geometry  model  used  in  CFD  simulation. 


temperatures  at  several  points  inside  the  reactor.  The  following 

assumptions  were  made: 

1.  Plug  flow  is  assumed  and  the  feed  and  oxidant  streams  are  pre¬ 
sumed  to  be  completely  mixed  at  the  reactor  inlet.  There  is 
neither  axial  nor  radial  dispersion.  This  hypothesis  of  complete 
mixing  was  verified  by  CFD  simulation  (Section 4.1). 

2.  The  model  assumes  that  the  oxidation  of  the  isopropyl-alcohol 
(IPA)  takes  place  in  a  single  step,  leading  to  the  formation  of 
carbon  dioxide  and  water  as  shown  in  Eq.  (1 ). 

C3H80  +  4.502  -*  3C02  +4H20  (1) 

3.  The  pressure  remains  constant  at  23  MPa  along  the  reactor. 

4.  Enthalpy  of  the  mixture  is  given  by  the  original  Peng-Robinson 
Equation  of  State  (PR-EoS)  [24].  Densities  are  calculated  by  the 
PR-EoS  with  constant  volume  translation  (VTPR-EoS)  [25].  In 
both  cases  van  der  Waals  mixing  rules  are  used.  The  volume 
correction  was  adjusted  for  each  pure  component  using  density 
data  from  the  NIST  database  [26]  and  from  Aspen  Plus  calcula¬ 
tions  using  the  SR-Polar-EoS  [27]  at  operating  pressure,  similar  to 
the  procedure  used  by  Lieball  [28].  Table  1  presents  the  volume 
correction  for  the  five  substances  that  compose  the  system.  Addi¬ 
tional  information  about  the  volume  correction  calculation  can 
be  found  in  Appendix  A.  The  translated  volume  of  the  mixture  is 
a  linear  function  of  the  component  mole  fractions. 


In  order  to  improve  the  optimization  process,  two  auxiliary  vari¬ 
ables  have  been  defined:  X\  =ln(/<0)  and  X2  =Ea/(R  •  850).  Also,  for 
convenience,  concentrations  were  substituted  by  mass  fraction. 
Thus,  Eq.  (5)  can  be  rewritten: 


r  =  exp 


X2  • 850  \ 
T  ) 


Wipa 


W02 

Mq2 


where  p  is  the  density  of  the  mixture.  The  value  850  is  a  charac¬ 
teristic  temperature  and  was  taken  in  order  to  give  X\  and  X2  the 
same  order  of  magnitude. 

The  optimization  procedure  uses  the  Matlab  function  lsqnon- 
lin,  which  consists  of  a  non-linear  least-square  solver  with 
Levenberg-Marquardt  algorithm  [30,31].  The  function  to  be  mini¬ 
mized  is: 


n 

fobj  =  -  ffc)2  (7) 

k= 1 

where  Tk  are  experimental  temperatures  and  Tk  the  calculated  tem¬ 
peratures.  These  temperatures  are  obtained  by  solving  the  set  of 
equations  for  wlPA  (Eq.  (2))  and  T  (Eq.  (4))  simultaneously  using 
Matlab  solver  ode 1 5s.  The  parameters  fitted  are  X\  and  X2  from 
which  the  Arrhenius  parameters  can  be  obtained. 


For  a  differential  control  volume  inside  a  plug  flow  reactor,  the 
material  balance  for  the  fuel  (IPA)  results  in  Eq.  (2): 

dwipA  _  r  •  MIPA  •  5 
dz  ~  m 

where  WiPA  is  mass  fraction  of  isopropyl-alcohol,  z  the  length  of  the 
tubular  reactor,  r  the  reaction  rate,  MIPA  the  molecular  weight  of 
isopropyl-alcohol,  S  the  cross  section  of  the  tubular  reactor,  and  m 
the  total  mass  flow. 

The  energy  balance  for  the  same  control  volume  leads  to  Eq.  (3). 
dh  AHr  r  S 

dz~  fn  1  J 

where  h  is  the  specific  enthalpy,  and  A Hr  is  the  heat  of  reaction, 
calculated  from  the  enthalpies  of  formation  of  each  component. 
Assuming  no  pressure  changes  dh  =  cpdT ,  Eq.  (3)  can  be  rewritten  as 
Eq.  (4). 


dT  _  A Hr  r  S 
dz  ~  m-Cp 


where  cp  is  the  specific  heat  capacity  of  the  mixture.  In  Eqs.  (2)-(4) 
the  reaction  rate  r  can  be  described  according  to  the  Arrhenius  law: 


r  -  k0 exp 


Reaction  order  regarding  IPA  is  assumed  as  equal  to  1,  as  sug¬ 
gested  in  previous  works  [29,12].  Oxygen  reaction  order  cannot  be 
assumed  as  0,  since  there  is  no  significant  excess  in  experimental 
data.  Therefore,  the  model  assumes  first  order  reaction  behavior 
for  both  fuel  and  oxidant.  The  assumed  order  of  reaction  gives  the 
S-shape  profiles  of  temperature  observed  in  experimental  data. 


3.2.  Study  of  the  mixing  process  of  water  and  air  flow 


The  assumption  of  instantaneous  mixing  taking  place  in  the 
tubular  reactor  was  verified  by  3D  simulation  of  the  tee-junction. 
The  extent  to  which  air  and  feed  flow  are  totally  mixed  was  inves¬ 
tigated  by  obtaining  the  concentration  profiles  in  different  lengths 
of  the  reactor.  In  these  simulations  a  non-reacting  flow  was  con¬ 
sidered.  The  geometry  of  the  model  includes  the  tee-junction  and 
a  25  cm  segment  of  the  reactor  and  the  metallic  wall,  as  shown  in 
Fig.  1.  The  domain  is  limited  by  the  longitudinal  plane  of  symmetry 
and  half  of  the  tee  and  tube  wall.  The  mesh  is  composed  of  436,000 
mixed  cells  at  the  fluid  region,  and  529,000  tetrahedral  cells  at  the 
solid  region.  The  fluid  region  mesh  is  formed  by  hexahedral  cells  at 
the  near  wall  region  and  tetrahedral  cells  at  core  region,  assuring 
any+  value  below  5. 

The  model  includes  the  five  species  that  constitute  the  system: 
water,  oxygen,  nitrogen,  carbon  dioxide  and  isopropyl-alcohol,  and 
solves  the  Reynolds-averaged  transport  equations  [33]  for  mass, 
momentum,  energy  and  species,  as  shown  in  Eqs.  (8)— (1 1 ): 


dp  3  ,  s.  n 

dt  dxi 

3  3 

3t  +  dx-  pUiUj  = 


dp  3  - 

dxi  dxj  1  J 


3  3  d  f  dT 

-(pE)+jrlui(pE  +  p)]  =  irl(k  +  Xt)ir  ]  +Sh 

dt  dXi  dxi  \  dxi 


3  3  3 

ol(PWk)  +  o-(pu/Wk)  =  - 
dt  dxi  dxi 


p(D  +  Df ) 


3  wk 

dxi 


+  R/, 


(8) 

(9) 

(10) 

(11) 
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where  is  the  fluid  velocity  in  the  direction,  p  the  pressure, 
E  is  the  total  energy  and  wk  the  mass  fraction  of  species  k.  The 
source  terms  Sh  and  Rk  represent  the  heat  of  reaction  and  the  net 
rate  of  species  k  production  by  chemical  reaction.  Thus,  for  the  non¬ 
reacting  flow  model,  their  values  are  zero.  The  viscous  stress  tensor 
r,j  is  shown  in  Eq.  (12): 


(12) 


where  pt  is  the  molecular  viscosity  calculated  as  a  mass-fraction 
average  of  the  properties  of  the  pure  components. 

In  Eq.  (10),  the  total  energy  E  is: 


r  ,  P  ud 

E  =  h  -  -  +  — 


P 


(13) 


Table  2 

Coefficients  of  polynomial  equation  fitted  for  enthalpy  of  reaction  of  IPA  oxidation 
as  function  of  temperature  at  23  MPa.  Temperature  in  °C,  enthalpy  in  J/mol. 


a0 

-2.0353  x  106 

b0 

-2.0780  x  106 

a\ 

-5.7254  x  10 

b\ 

9.7755  xlO2 

a2 

4.5344 

b2 

-1.7952 

a3 

-2.2537  x  10~2 

b3 

1.4322  xlO”3 

a4 

3.4218  xl0“5 

b4 

-4.1826  x  IQ”7 

reaction  chamber  and  transpiring  chamber.  A  complete  description 
of  the  base  model  is  given  by  Bermejo  et  al.  [23]. 


4.  Results 

4.2.  Study  of  mixing  influence 


In  Eqs.  (8)-(13)  p  is  the  mixture  density  calculated  by  VTPR- 
EoS,  X  the  thermal  conductivity  calculated  by  TRAPP  method  for 
high  pressure  fluids  [33],  D  the  diffusion  coefficient  calculated  as 
the  mass-fraction  average  of  the  pure  species,  and  specific  enthalpy 
h  is  given  by  PR-EoS.  The  turbulent  thermal  conductivity  Xt,  the 
turbulent  diffusion  coefficient  Dt  and  the  Reynolds  stress  tensor 
-pu'm'j  are  defined  according  to  the  turbulence  model  being  used. 
In  this  work,  the  Realizable  k-e  model  has  been  used  in  order  to 
model  turbulence,  as  suggested  by  Sierra-Pallares  et  al.  [34],  with 
enhanced  wall  treatment. 

The  model  was  solved  using  the  commercial  CFD  software 
ANSYS  FLUENT  1 2.0,  with  segregated  pressure  based  solver,  SIMPLE 
algorithm  for  pressure-velocity  coupling,  and  second  order  upwind 
discretization  for  energy,  momentum,  k-e  and  species  equations. 
The  boundary  conditions  for  the  fluid  region  are  mass  flow  inlet, 
pressure  outlet,  symmetry  and  no-slip  condition  (at  the  interface 
with  the  solid  domain). 

The  test  case  was  modeled  also  with  other  near  wall  treatments, 
like  the  low-Reynolds  k-e  turbulence  model  [35],  and  meshes  with 
different  refinement  levels,  but  showed  no  significant  differences 
respect  the  model  described  here. 

3.3.  Transpiring  wall  reactor 


Fig.  2  shows  the  dispersion  of  species  according  the  result  of  a 
non-reacting  simulation  of  the  tee-junction.  Concentrations  (in  ver¬ 
tical  axis)  are  normalized  by  the  mass  fraction  of  equilibrium,  when 
the  species  are  completely  mixed.  The  simulation  corresponds  to 
an  experiment  using  a  feed  flow  of  5.9  kg/h  with  4%  mass  isopropyl 
alcohol  at  488  °C  and  3.7  kg/h  of  air  at  389  °C,  both  at  23  MPa,  which 
gives  a  Reynolds  number  of  26,000. 

For  these  conditions  the  simulation  shows  that  IPA  concentra¬ 
tions  is  stabilized  nearly  10  cm  after  the  mixing  point.  In  our  facility, 
the  interval  between  two  temperature  measuring  points  was  about 
20  cm  and  the  total  length  of  the  reactor  was  1 .5  m,  making  insignif¬ 
icant  the  mixing  length.  Thus,  the  hypothesis  of  instantaneous 
mixing  is  acceptable  and  the  reactor  can  be  modeled  as  plug  flow. 
Also,  a  visual  reference  of  the  mixing  process  can  be  seen  on  detail 
in  Fig.  2,  where  the  simulated  contours  of  IPA  concentration  are 
presented. 

4.2.  Kinetic  model 

The  fitted  parameters  and  their  confidence  interval  (Cl)  are  pre¬ 
sented  in  Eq.  (15).  The  confidence  interval  of  the  parameters  was 
calculated  by  the  Matlab  function  nlparci  which  calculates  the 
95%  confidence  intervals  for  parameters  in  nonlinear  regression. 


A  model  previously  developed  for  describing  the  behavior  of 
the  transpiring  wall  reactor  [23]  was  modified  in  order  to  describe 
the  performance  of  the  reactor  with  a  hydrothermal  flame  inside, 
incorporating  the  kinetic  model  obtained  in  this  work.  The  model 
assumes  stationary  conditions  considering  only  mass  and  energy 
balances,  resolved  by  the  Runge-Kutta  method  of  4th  order,  while 
the  momentum  equation  was  neglected.  The  pressure  along  the 
reactor  was  considered  constant  and  equal  to  23  MPa.  Enthalpy 
and  heat  capacities  were  calculated  using  PR-EoS.  Densities  were 
calculated  by  VTPR-EoS  considering  the  exact  composition  of  the 
mixture.  Global  heat  transmission  coefficients  were  calculated  as 
a  function  of  the  total  mass  flow,  and  the  coefficient  for  the  heat 
transmission  of  the  outer  wall  with  the  environment  was  taken 
as  0.5Wm_2I<_1.  The  heat  of  reaction  for  IPA  oxidation  were 
calculated  for  different  temperatures  using  a  polynomial  fit  (Eq. 
(14)). 

f  a0  +  a{Y  +  a2T2  +  a3T3  +  a4T4,  T  <  404 °C 
AHr(T)  =  <  (14) 

[  b0  +  b{Y  +  b2T2  +  b3T3  +  b4T 4,  T  >  404 °C 

where  the  coefficients  of  the  polynomials  are  presented  in 
Table  2. 

This  model  assumes  the  injector  as  a  PFR,  the  upper  part  of  the 
reaction  chamber  as  a  CSTR,  and  the  reaction  chamber  is  modeled 
as  a  plug  flow.  The  model  includes  heat  exchange  between  injec¬ 
tor  and  reaction  chamber,  and  heat  and  mass  exchange  between 


Reactor  length  (m) 

Fig.  2.  Maximum  and  minimum  species  concentration  relative  to  equilibrium  con¬ 
centration  in  the  tubular  reactor,  showing  the  dispersion  along  reactor  length:  IPA 
(black  line),  O2  (grey  line).  The  central  dotted  line  marks  the  equilibrium  value. 
Above,  detail  of  simulation:  IPA  (black),  02  (white). 
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Reactor  length  (m) 

(b) 


Reactor  length  (m) 

Fig.  3.  Comparison  of  experimental  temperatures  (symbols)  along  the  tubular  reac¬ 
tor  and  model  predictions  (continuous  line)  using  the  fitted  kinetics.  Dashed  lines 
indicate  the  predictions  in  the  limits  of  Cl.  (a)  5.9  kg/h,  4%IPA  at  395  °C,  air  as  oxidant; 
(b)  1 0.6  kg/h,  2%IPA  at  404  °C,  O2  as  oxidant. 


Thus,  the  kinetic  parameters  are  shown  in  Eq.  (16): 

k0  =  (9.308  ±  3.989)  x  107(m3  s”1  kmol-1 ) 

,  (16) 

Ea  =  89.441  ±  2.457(kJ  mol~ 1 ) 

Fig.  3  compares  the  experimental  temperature  profile  with  the  pre¬ 
diction  of  the  model  using  the  parameters  found  in  Eq.  (16).  The  data 
corresponds  to  experiments  at  two  different  inlet  temperatures, 
flows  and  type  of  oxidant  used:  air  or  oxygen. 

The  model  prediction  agrees  with  experimental  data,  as  the 
experimental  temperatures  are  inside  the  model’s  interval  of  con¬ 
fidence.  The  least  square  error  over  all  experimental  points  used  to 
perform  the  fitting  is  10.8%. 

4.3.  Parametric  analysis 

The  kinetic  model  was  applied  in  a  parametric  analysis  to  study 
the  influence  of  inlet  temperature  and  fuel  concentration  on  flame 
formation.  Once  again,  a  plug  flow  reactor  was  assumed  and  the 
results  are  showed  as  function  of  residence  time. 

Fig.  4  shows  how  the  fuel  concentration  affects  the  flame 
formation.  Three  different  concentrations  where  tested  at  inlet 
temperature  of  350  °C,  and  air  in  10%  excess.  The  analysis  of  the 
reaction  rate  shown  in  Fig.  4(a)  indicates  that  a  feed  with  3%  IPA  at 
350 °C  is  slowly  oxidized,  while  a  feed  with  5%  IPA  is  rapidly  oxi¬ 
dized  inside  a  narrow  zone  that  could  be  recognized  as  a  flame. 
Fig.  4(b)  shows  the  temperature  rising  for  the  same  cases  and  sug¬ 
gests  a  flameless  oxidation  for  3%  IPA  and  a  flame  oxidation  for  5% 
IPA.  The  concentration  profiles  are  shown  in  Fig.  4(c)  where  can 
be  seen  a  clear  separation  of  reagents  and  products  for  the  5%  IPA 
feed.  These  results  are  in  accordance  with  the  work  of  Serikawa 
[21],  that  found  stable  hydrothermal  flames  using  a  tubular  reac¬ 
tor  for  IPA  oxidation,  when  IPA  concentrations  were  above  4%  (or 
5  vol%,  since  volume  fraction  was  used  in  that  work  instead  of  mass 
fraction). 

The  effects  of  inlet  temperature  are  shown  in  Fig.  5.  Feeds  with 
5%  IPA  and  air  with  1 0%  excess  at  three  different  injection  tempera¬ 
tures  are  compared.  Fig.  5(a)  shows  that  the  injection  temperatures 


(b) 


Fig.  4.  (a)  Reaction  rate,  (b)  temperature  and  (c)  species  mass  fractions  as  function 
of  residence  time  in  plug  flow  reactor  at  injection  temperature  of  350  °C  for  different 
inlet  IPA  concentrations:  3%  (continuous  line)),  4%  (dotted  line)  and  5%  (dashed  line). 


Residence  time  (s) 

Fig.  5.  (a)  Reaction  rate,  (b)  temperature  as  function  of  residence  time  in  plug  flow 
reactor  at  5%  IPA  for  different  temperatures  of  injection:  330  °C  (continuous  line)), 
350 °C  (dotted  line)  and  370 °C  (dashed  line).  Grey  line  corresponds  to  ignition 
temperature. 
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have  influence  on  flame  position.  Reducing  injection  temperature 
from  370  to  330  °C  causes  a  delay  of  approximately  0.2  s  in  the  flame 
ignition.  With  the  characteristics  of  our  experimental  plant,  this 
delay  is  equivalent  to  a  displacement  of  60  cm.  Fig.  5(b)  shows  tem¬ 
perature  profiles  for  the  three  cases  and  the  ignition  temperature 
found  by  Serikawa  [21].  This  temperature  (470  °C)  corresponds  to 
the  beginning  of  the  maximum  slope  region  at  temperature  profiles 
predicted  in  this  work. 

4.4.  Simulation  ofTWR 

The  result  of  the  transpiring  wall  reactor  simulation  shown  here 
corresponds  to  a  piece  of  experimental  data  using  a  feed  flow  of 
27.7  kg/h  with  8%  mass  isopropyl  alcohol  and  23.9  kg/h  of  air,  an 
injection  temperature  of  246  °C  and  a  pressure  of  23  MPa.  The  tran¬ 
spiring  flow  was  16.9  kg/h  of  water  at  25  °C. 

The  results  of  this  simulation  show  that  the  model  correctly 
predicts  the  heating  and  flameless  oxidation  inside  the  injector 
and  also  the  flame  formation  at  upper  part  of  reaction  chamber.  In 
presence  of  a  hydrothermal  flame,  complete  oxidation  occurs  when 
injection  temperatures  are  as  low  as  246  °C.  Inside  the  injector,  the 
reagent  stream  is  heated  up  to  382  °C  due  heat  transmission  from 
the  reactor  chamber  through  the  injector’s  wall,  and  slow  oxidation. 
Outside  the  injector,  the  reactive  mixture  reaches  a  temperature  as 
high  as  709  °C  (experimental  value  was  694  °C).  The  temperature 
profile  predicted  by  the  model  is  compared  to  the  experimental 
temperature  data  measured  inside  the  injector  and  at  the  reaction 
chamber  in  Fig.  6(a).  A  good  agreement  with  the  experimental  data 
can  be  observed  here. 

Fig.  6(b)  shows  the  same  experiment  reproduced  by  the  pre¬ 
vious  model  of  our  group  [23].  The  kinetic  model  adopted  in  that 
previous  work  was  able  to  reproduce  experimental  data  for  inlet 
temperatures  of 300  °C  and  above,  but  failed  at  lower  injection  tem¬ 
peratures.  Other  models  existing  in  literature,  fitted  for  systems  at 
lower  temperatures,  are  not  fast  enough  to  correctly  predict  the 
rapid  reaction  rate  at  flame  zone.  Nevertheless,  the  model  devel¬ 
oped  in  this  work  shows  satisfactory  results  for  a  wide  range  of  inlet 
temperatures  and  concentrations. 


Temperature  (°C) 

Fig.  A.7.  Density  fitting  for  pure  species  using  PR-EoS  at  230  bar. 

5.  Conclusions 

A  new  simple  kinetic  model  for  supercritical  water  oxidation  of 
isopropyl-alcohol  in  hydrothermal  flame  regime  was  determined 
by  adjusting  the  temperature  profile  in  tubular  SCWO  reactors.  The 
model  is  able  to  predict  temperature  profiles  in  tubular  reactors  at 
different  IPA  concentrations,  feed  flows  and  using  air  or  oxygen  as 
the  oxidant. 

A  3D  CFD  model  confirmed  the  hypothesis  that  the  mixing  pro¬ 
cess  occurs  in  the  first  centimeters  of  the  tubular  reactor.  Thus,  the 
initiation  of  the  reaction  is  not  affected  by  the  mixing  process,  being 
the  reaction  controlled  by  the  kinetics. 

Given  the  inlet  conditions,  the  model  can  estimate  if  the  oxida¬ 
tion  will  be  flame  or  flameless  regime,  and  is  capable  of  reproducing 
both. 

This  kinetic  model  was  applied  to  a  model  for  the  transpir¬ 
ing  wall  reactor  developed  by  the  University  of  Valladolid  and  the 
experimental  data  were  successfully  described.  According  to  the 
results,  in  the  presence  of  a  hydrothermal  flame,  the  reactor  can  be 
maintained  in  steady  state  regime,  with  total  destruction  of  fuel, 
even  for  subcritical  feed  injections  (246  °C). 

The  new  kinetic  model,  allied  with  modeling  and  simulation, 
allows  for  a  better  understanding  of  SCWO  reactors  and  can  lead  to 
improving  their  design  and  optimization. 
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Fig.  6.  Experimental  (symbols)  and  model  (lines)  temperature  profiles  ofTWR  sim¬ 
ulation:  (a)  this  work,  (b)  previous  model.  Black  color  corresponds  to  injector;  grey 
color  represents  reaction  chamber. 
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Appendix  A.  Volume  translation  fitting 

Fig.  A.7  shows  the  fitting  of  the  densities  of  pure  compounds 
by  VTPR-EoS.  The  continuous  line  corresponds  to  NIST  data  as  a 
function  of  temperature  at  230  bar,  and  the  dotted  line  corresponds 
to  the  predicted  densities.  The  volume  correction  was  estimated 
for  each  substance  in  order  to  minimize  the  differences  between 
the  reference  data  and  the  predicted  values.  For  isopropyl-alcohol, 
density  data  are  not  available  on  the  NIST  database,  so  the  reference 
data  were  obtained  by  SR-Polar-EoS  calculations  in  Aspen  Plus. 
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